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Abstract: We consider the one dimensional totally asymmetric simple exclusion 
process with initial product distribution with densities < po < p\ < ... < p n < 1 
in (—00, ci£ _1 ), [ci£ _1 , C2£ _1 ), . • • , [c n £ _1 , +00), respectively. The initial distribu- 
tion has shocks (discontinuities) at e~ 1 Ck, k = l,...,n and we assume that in 
the corresponding macroscopic Burgers equation the n shocks meet in r* at time 
t*. The microscopic position of the shocks is represented by second class particles 
whose distribution in the scale e -1 / 2 is shown to converge to a function of n inde- 
pendent Gaussian random variables representing the fluctuations of these particles 
"just before the meeting". We show that the density field at time in the 

scale £ -1 / 2 and as seen from e -1 r* converges weakly to a random measure with 
piecewise constant density as e — > 0; the points of discontinuity depend on these 
limiting Gaussian variables. As a corollary we show that, as e — > 0, the distri- 
bution of the process at site e~ l r* + e~ 1 / 2 a at time e~ l t* tends to a non trivial 
convex combination of the product measures with densities pk, the weights of the 
combination being explicitly computable. 
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§1. Introduction and results. 

It is well known that the hydro dynamical behavior of the one dimensional asym- 
metric exclusion process is described by the inviscid Burgers equation 



where 7 is the mean of the jump distribution. Since equation (1.1) develops discon- 
tinuities, one has to be careful about the precise statement, but loosely speaking, if 
(r, t) is a continuity point of p(r, t), for a given initial measurable profile po( • )> then 
at the macroscopic point (r, t) the system is distributed according to the measure 
v p{r,t)i where v p is the product Bernoulli measure on {0, 1} Z with v p (r}(x) = 1) = p 
for all x. This is known as local equilibrium. As it is known, the exact statement 
involves a space/time change, under Euler scale, and for all these developments we 
refer to Andjel and Vares (1987), Rezakhanlou (1990), Landim (1992). 

The problem with which we are concerned here involves the description of the 
(microscopic) behavior of the system at certain discontinuity points (r, t) (or shock 
fronts) of the solution of equation (1.1). For example, if 7 > and the initial 
profile is a step function po(r) := al{r < 0} + (3l{r > 0}, with < a < (3 < 1, 
the entropy solution of equation (1.1) is p(r, t) = al{r < vt} + (3\{r > vt}, where 
v := 7(1 — a — (3) is the velocity of the shock front and 1{-} is the indicator function 
of the set {•}. This description is valid only for continuity points. The investigation 
of what happens to the system if looked from this shock front was first studied by 
Wick (1985) for a different model, and for the asymmetric simple exclusion in the 
particular situations a = and a + (3 = 1 by De Masi et al (1988) and Andjel, 
Bramson and Liggett (1988), respectively. They all proved that at the shock front 
one sees a fair mixture of v a and v$. This result has then been extended so as to 
cover all cases 0<a<,9<lby Ferrari and Fontes (1994), from now on referred as 
[FF]. [FF] worked with the nearest neighbor asymmetric exclusion process, whose 
generator is the closure of 



(1.1) 



d t p + 1 d r (p(l-p)) = 



(1.2) 



xEZy=x±l 
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for / a cylinder function in {0, 1} Z , with 




i](z) iiz^x,y 
rj(y) if z = x 
77(a) i£z = y, 



where p(x,x + 1) := p, p(x,x — 1) := q := 1 — p, with 1/2 < p < 1. This process 
was first studied by Spitzer (1970). Calling /ia^ the product measure on {0, 1} Z 
with site marginals 



denoting S t as the semigroup corresponding to the above generator and 6 X as the 
space shift /J.0 x (f) := / f{9 x rf)ii{drf), with 6 x r](z) := rj(x + z), [FF] proved that 



as t — > +00, and where [x] denotes the integer part of x. This corresponds to 
the exact statement made above, under Euler scale, and for the macroscopic point 
(r, t) in the front line, i.e. r = vt. In fact, in the above mentioned references, more 
detailed analysis is performed, by looking at the microscopic structure of the shock 
represented by a second class particle. A second class particle jumps to empty 
sites with the same rates as the other particles, but interchanges positions with the 
regular particles at the rate holes do. A formal definition using coupling is given 
in the next section. Calling X t the position of a second class particle added at 
the origin, the process as seen from the second class particle 6x t Vt has distribution 
asymptotically product to the right and left of the origin with densities a and (3 
respectively, uniformly in time. The velocity of the second class particle is the same 
as the velocity of the shock in the Burgers equation: E fla Xt = 7(1 — a — f3)t. [FF] 
proved that the fluctuations of the position of the particle are Gaussian: calling 



(1.3) 




(1.4) 



2^ + ^) 



X t := X t - vt, 



(1.5) 



Xt/Vi 
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where Af a ,/3 is a centered Gaussian random variable with variance r y((3 — a)~ 1 (a(l — 
a) + (3(1 — f3)). With this result in hand [FF] proved that, if — oo < a < +oo, the 

distribution of the process at time t at the point vt + ay/i converges to a mixture 
of v a and vp\ more precisely, for real a, 

(1.6) Ha,/3StO[ vt +ay/i\ VaP{Na,l3 > a) + VpP(jV a ,p < a) 

which in particular yields \{y a + vp) for a = and interpolates between v a and vp, 

as a varies from — oo to +oo. 

Our goal is the consideration of two or more shock fronts and the description 

of the microscopic behavior of the system at the (macroscopic) time and position 
of their meeting. To avoid unnecessary technical difficulties, we look the totally 
asymmetric case: we assume 

p = 1, q = which implies 7 = 1 in equation (1.1). 

The extension to ^ < p < 1 will be briefly discussed at the end. 
We consider points Cfc, densities pk and the existence of a space-time point (r*, £*) 
such that 

(1.7 ) - oo = c < ci < • • • < c n < c n+1 = oo, 

(1-8) 0<p O <---<Pn<l, 

(1.9) r* = c k + (1 - p k -! - p k )t*, k=l,...,n. 

With this assumption the entropy solution to equation (1.1) with initial data 

n 

(1.10 ) p(r) := Pk !{cfc < r < c k+1 )}, 

k=0 

has the property that all n shocks meet at r* at time t* . More precisely, the entropy 
solution is given by 



(1.11) 



n 

p(r,t) = ^2p k l{c k (t) < r < c k +i(t))} 

k=0 
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where 

(1.12) c k (t) = 



c fc + (1 - p k -i - Pk)t for t < t* 

r* + {1- po- p n )(t-t*) foTt>t* 



Notice that after t* only the extreme densities po and p n are seen. 

In our discussion we shall need to consider entropy solutions starting with more 
general increasing step profiles, i.e. not necessarily with all shocks meeting at 
(t*, r*). For this let (p k ) satisfy (1.8), — oo = bo < b± < ... < b n < b n+ i = +oo and 
consider the entropy solution to the Burgers equation with initial data 

n 

(1.13 ) A(r) := Pk Hh < r < b k+1 )}, 

k=0 

which is given by 

n 

(1.14) X(r,t) = J2p k l{b k (t)<r<b k+1 (t))} 

(1.15) = l-X+(b k (t),t)-X-(b k (t),t) 

6fc(0) = ft*,, i = l,...,n. 

where A^(r, t) are the right and left limits of A(r, t): A^(r, t) = lim r /^o A(r ± r', t), 
with r' > 0. In words, b±(t) < ■ ■ ■ < b n (t) represent the shock fronts, which move 
initially as bj(t) = bj+t(l— Pj-i~ Pj) for 1 < j < n, until two or more of them meet. 
When this happens the involved shock fronts coalesce, with the disappearance of 
all intermediate densities and the front keeps moving with a new velocity given by 
one minus the sum of the two densities, to its left and to its right, i.e. the two 
densities which form the shock. This simple description, mainly due to the fact 
that we are in the one dimensional situation and the initial profile is an increasing 
step function, allows to define the following map ip : R n — > M n , ip = (ip k ), which 
will be fundamental in the determination of the distribution of the process at the 
macroscopic meeting point of the shocks: 

Definition of if;. Given a vector x = (x±, . . . , x n ) in M n let us take a time t(x) 
large enough so that defining 

(1.16 ) b k (x) := x k - t(x) (1 - p k - X - p k ) 
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then 

(1.17) h{x) <■■■ <b n {x). 

That is, if we consider a family of one shock solutions A fc (r, t), 1 < k < n starting 
with 

X k (r) := p k -il{r < b k (x)} + p k l{r > b k (x)} 

then x k is the position of the shock at time t(x), for the k-th equation. What we 
denote by ip(x) is then the position of the shock fronts at time t(x) in equations 
(1.14)-(1.15) starting with (1.13) for (bk = bk(x)). It is easy to see that the definition 
is well posed, i.e. it does not depend on the value of t(x) provided (1.16) and (1.17) 
hold. 

The coordinates of the vector ip(x) = (V'fe(^)) are convex combinations of some 
Xj. In the particular case of n = 2 this becomes 

{x k if x\ < x 2 

P2 - PO Pi - PO 

for k = 1,2. 

We consider a family of product measures p e on {0, 1} Z with profile p: the 
marginal one-site distribution is given by 

p e (r](x) = 1) = p(ex), 

where p(-) is given in (1.10) and corresponds to the case when all shocks meet at 
the same point r* at time t*. Considering pf as the initial measure, our first goal 
is to look at the asymptotic distribution of the process at time r e around site [x e ] 
in the scale £ -1 / 2 as e — * 0, where 

(1.19) T e :=t*e-\ x s :=r*e-\ 

Let 3^o := — oo, y n +i '■= oo and 



(1.20) 



(yi,...,y n ) :=^{Xi,...,x n ) 
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where (X\, . . . , X n ) are independent centered Gaussian random variables with vari- 
ances (Di, . . . , D n ) given by 



^ 21 ^ Dk . = Pfc-i(l -Pfc-i) + Pfc(l -Pfc) 

Pk-Pk-i 



t* 



Theorem 1.1. Let jU e (0 < e < 1) be the product measure on {0, associated 
to the profile p(-), given by (1.10). Let St denote the semigroup associated to the 
totally asymmetric n.n. simple exclusion process, corresponding to the generator 
given by (1-2) with p = 1. Then, for a£l, 

n 

(1.22) f* e £ t . e -i0 [r . e -i+«-i/*] £^ X> Pfc P(y*<a<y fc+1 ) 

~* fc=0 

where 3^fc are defined in (1.20). 

We shall see that, as in (1.5), (3^i, y n ) represents the limiting fluctuations of 
the shocks at time around r*£ _1 . In the case of only one shock (n = 1) 

this agrees with (1.6). The crucial difference is that if, say, the k-th and (k + l)-th 

microscopic shocks meet before they coalesce and change the velocity and 

the intermediate zone of density pk disappears. This explains the weight of ^ Pfc in 

(1.22) : it is the same as the probability that (a) the k-th and (k + l)-th shocks 
have not collided yet and (b) these shocks are to the left and right of the point 
we are looking at, respectively. The coalescing dynamics of the microscopic shocks 
in the scale £ -1 / 2 is as in the Burgers equation. Its relation with the n one-shock 
dynamics — represented by (Xk) — is given by the function tfj. 

We turn now to the profile seen from the meeting point and time of the shock 
fronts, scaled by e~ 1 / 2 . For e > and a local function /, consider the (random) 
measures on the real line given by 

(1.23) A £ (da) := e^Ysfi^+l^VrJ^Ada), 

xei, 

where x e and r e are given by (1.19), 5 e i/2 x is the Dirac delta measure at e 1 l 2 x, and 
rj t is the configuration at time t for the initial measure . 
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Let (3^i, • • • , y n ) be as in (1.20) and consider the random measure 

n 

(1.24) A(da) :=J2" Pk (f) 1{34 < a < y k+1 } da 

k=0 



Theorem 1.2. A g converges in law to A as e — > 0, for the usuai weak topology on 
the space of measures. 

The analysis is based on the well known strategy of identifying microscopically 
the shocks with second class particles, defined through the so called basic coupling 
of different versions of the process starting with measures with different uniform 
densities. From the dependence of their locations on the initial condition, which is a 
one-shock fact as in [FF], we can ascertain their distribution around the macroscopic 

meeting place [x e ] at the macroscopic meeting time r e ; in the scale e -1 / 2 this is 
given by (X\, . . . , X n ). If we look at the positions of the shocks at time r e — ae -1 / 2 , 
the law of large numbers apply. This is also a one-shock fact (Ferrari (1992)). If 
a is big enough these positions are ordered, with very large probability, and they 
are represented by the variables bk(Xi, . . . , X n ) as in (1.16), taking t(x) = ae~ 1 / 2 . 
From then on, using the invariance of the product measures involved, we can follow 
their trajectory up to the meeting time through successive applications of the one- 
shock law of large numbers. The final positions of the shocks in the scale e~ 1 / 2 are 
given by the variables 34 + [x e ]. 

Theorem 1.2 follows from the weak convergence of a suitable function of the 
second class particles to the variables 34, and the asymptotic properties of the 
measure as seen from the second class particles. We show Theorem 1.1 as a corollary 
of Theorem 1.2, by using the attractiveness of the process. The proof avoids proving 
firstly the translation invariance of the weak limits of Theorem 1.1, as in the proof 
given by [FF] for the one-shock case. 

Acknowledgment: The authors are indebted to an anonymous referee for 
bringing in the discussion of Theorem 1.2 and for noticing that our first proof 
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of Theorem 1.1 was indeed yielding this profile description. 
§2. Second class particles and more. 

Ferrari, Kipnis and Saada (1991) and Ferrari (1992) have shown that in the case 
that /j,a : p is the initial measure, the shock front is well described by X t , the position 
of a "second class particle" initially located at the origin. This fact together with 
the validity of a central limit theorem for X t as t — > +oo, proven by [FF], are 
the essentials for the proof that in the totally asymmetric case fi> a ,pSt0[(i-a-i3)t] 
tends to \(v a + vp). For all this, as well as in the present work, the main tool is 
coupling. To realize a coupling of several evolutions of the simple exclusion process, 
corresponding to several initial configurations, is particularly simple through the 
graphical construction: to each pair of sites (x,x + 1), let us associate a Poisson 
process with rate 1, and at each of its occurrences we put an arrow x — > x + 1. 
Construct all such Poisson processes as independent in some space (ft, A, P). Given 
any realization of arrows u> and a initial configuration r\, we may realize an evolution 
r\t corresponding to L imposing that whenever an arrow x — > x + 1 appears, if there 
is a particle at x and no particle at x+1, then this particle moves to x+1; otherwise, 
nothing happens. (In the non totally asymmetric cases, we have arrows (x, x + 1) 
with rate p and (x, x — 1) with rate q, all Poisson processes being independent.) 
Since the probability of two simultaneous arrows is zero this construction makes 
sense and defines the process of interest. (One could also state the coupling via a 
suitable generator, cf. Liggett (1976) and Chapter VIII of Liggett (1985).) From 
this coupling, the attractiveness property of the dynamics is immediate: if rf and 
rj 1 are two initial configurations such that rf < r] 1 (i.e. rf(x) < ^(x) Wx) and we 
write 77°, r)\ for their evolutions using the same arrows, then rfl < r)\ for each t. 

Let us then couple in this way realizations of the exclusion process with random 
initial configurations, distributed according to v Pk , k = 0, . . . , n. In fact we have 
a small perturbation since we will add particles at [£ _1 Cfc]. For example, let U = 
(U x )xei De a family of i.i.d. random variables uniformly distributed on [0, 1], taken 
as independent of all the Poisson processes of arrows, where, if needed, we enlarge 
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the basic probability space (Cl,A, P). Then, for (c/~) and (pk) as in Section 1 define 



Using the same graphical construction above described we may consider the 
simultaneous evolution of all these configurations, which we denote by a f fe , k = 
0, ...,n on the space {£l,A,P). The marginal distribution of a\ is the simple 
exclusion process under the invariant distribution v Pk . 

Consider the configurations on {0, 1} Z given by $ = of and for k = 1, . . . , n, 



It is easy to see that when considering the joint motion of . . . , then for 
j < k the £ J particles have priority over the £ fc particles: if there is a £ J particle 
at site x, a £ fc particle at x + 1 and an arrow from x to x + 1, then they exchange 
positions. Otherwise, the interaction is the usual exclusion. 

Denote X k the position at time t of the £ h particle which was at site [c^e -1 ] at 
time 0. 

The essential tools in [FF] (with as the initial measure) include the joint 
realization of the evolutions (rj' t , Z t ) and (cx° , X£), where r/' is distributed as \x a ^ 
conditioned to have the site occupied by a second class particle, and Z t describes 
the position of this single second class particle, while oy, and X\ were defined 
above for n = 1, po = a, p\ = (3. 

When p = 1 the above coupling can be achieved by 



and Z t = X}. Together with the law of large numbers and the central limit theorem 
for X}, this coupling is the basic ingredient in [FF]. Still restricting ourselves to 
the case p = 1, a natural extension of this to our evolution rjt starting with p e is 
the following. 



<T k ([e- 1 c j ]) = l{j<k}, j = l,...,n, 

a k (x) = 1{U X < p k }, xeZ\ {[e-^i], . . . , [e' 1 ^]}. 



(2.1) 



tf(x)=a*(x)-at 1 (x) 
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Definition of (Y t k ). Let us recall that we consider the case p = 1, and X k are as 
just described above. Let F t ° = — oo, Y™ +1 = oo. For fe = 1, . . . , n, we define y/ 1 ' 1 
and ti inductively in % > as follows. Let to = and Y t k '° = X k for all k = 1, n 
and t > 0. Having defined tg and f/ '^ for all k = 1, n and t > t^, if ^ = oo, we 
stop the inductive procedure; otherwise, let tg + i = inf{t > tg : Y^' £ = Y^ +1 ' £ + 1} 
(with the usual convention that inf0 = oo). Denote by ig the index involved. If 
finite, tg+i is the time of the first crossing after time tg of two particles whose 

positions at time tg are Y t f. At time tg + i, a a lt ~ l \a lt+1 discrepancy appears at 
the position Y^^ 1 . Since p = 1, from tg + i on, Y^ e ' e and y t ^ +1 ^ do not uncross, 
so we can ignore the particles and consider only the evolution of (a 1 : i ^ ig). 
Accordingly, we define, for t > tg + i, Y^ e ' £+1 = Y* e+1,£+1 to be the position of this 
a ii-i\ a t e +i discrepancy at time t. We also make Y^' £+1 = Y^' 1 for i ^ ig,ig + 1. 
This concludes the inductive step. 

Notice that the inductive steps in the above procedure correspond to instants of 
crossings of shocks, after each of which we coalesce the crossing shocks into a single 
one, which is then made to follow a new discrepancy, leaving the remaining shocks 
untouched. 

We are ready to define (Y t k ). For all 1 < k < n and t > 0, let Y t k = Y t k/ , if 
tg <t < tg + i. In case t n < oo, t n+ i is defined as oo. 
We further define 

(2.2) rj' t (x) =j^a k t {x)l{Y t k < x < Y t k+1 }. 

fc=0 

Notice that some of the indicator functions may vanish. This happens when the 
corresponding Y particles coalesce before t. Since p = 1, rf t has the same law as the 
evolution starting from \i e conditioned to have the sites [c^e -1 ] occupied by second 
class particles with respect to the other particles of ?]' such that the second class 
particles of lower labels have priority over the ones with higher labels. 



Theorem 2.1. Writing = X* e - [x £ ] and Y^ = Y^ - [x £ ], where r e and x £ 
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are given by (1.19), we have: 

(2.3) lime 1/2 (J^ ,,...,*») = . . . , * B ) 

(2.4) hW/ 2 ^,...,^) £ (^,...,^) 

where (Xk) are i.i.d. centered Gaussian with variances D k given by (1.21) and 
(3^i, • • • , 3^n) = ip(Xi, . . . , A^i), with •0 defined in Section 1. 

Before proving Theorem 2.1, let us recall that from Theorem 1.1 of [FF] we have 
(2-5) e^Xl ^X k , 

with (X*) and (X k ) as in the previous statement. We discuss independence below. 
The basic idea to show this theorem is to look at the system "just before" r e (in 
macroscopic time); more precisely: let a > and = r e — ae -1 / 2 . We shall make 
e — > and then a — > +oo. 

Presumably, if e is small, but a is large enough, with overwhelming probability 

the second class particles in the n one-shock systems, (X^ a , ■ ■ ■ , X™ a ) have not 
yet crossed each other and so FA = Xtl a for k = l,...,n. This is due to the 

£ £ 

control on the asymptotic behavior of Xj°, known from [FF] and a reasoning as 
the one used to define b k {x) in (1.17). This allows to compare the system with n 
independent systems of only one shock each, and we conclude that (Y^ a , . . . , Y™ a ) 
are at asymptotically independent Gaussian distances from their expected values 
(on the £ -1 / 2 scale). (Notice that since p = 1, then each pair X 1 and X^ with i < j 
cross each other only once.) 

Since, furthermore, (Y* a , . . . , Y™ a ) are at distance of order e -1 / 2 apart from each 

other, in the time interval from to r £ (of length ae -1 / 2 ) the fluctuations are 
negligible (we are in the regime of the hyperbolic scaling, thus of the law of large 
numbers for Y k ) and it is as if we have a completely deterministic situation. 

Rigorous versions of the facts described on the previous two paragraphs lead 
straightforwardly to Theorem 2.1. We state and prove them in the paragraphs 
below. 
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The first step is to describe the behavior of (Xl a , . . . , X%). Since each of the 

£ £ 

coordinates X k correspond to the motion of a second class tagged particle in a 
single shock situation, the essential is their dependence on the initial condition. 
Indeed, if 

l(c k +t(pk-Pk-i))/e] [c fe /e] 

x=[c k /e] + l x=[(c k -t(p k -p k _ 1 ))/e] + l 



it follows from [FF] that 



(2.6) 



et^E 



vk 
A t/e 



Ck 



Pk — Pk-i 



Remark. Indeed the result in [FF], for X k , is for uniform product measures off 
the initial location of the tagged particle. A straightforward conditioning argument 
shows it to be true for any finite perturbation of that measure. From Theorem 3.1 
in (Ferrari (1992)), but also from (2.6) it holds 

(2.7) lim eX* lf = c k + (1 - p k -i - p k )t, in probability. 

Proof of (2.3). If t < t* the random variables N k e are independent since by (1.9) 
they are functions of the initial configuration at disjoint sets of sites. This and the 
product nature of the initial distribution allows us to use the central limit theorem 
for sums of Bernoulli random variables to show (2.3). □ 

To compute the limiting joint probabilities of (2.4) let A > and consider the 
following intervals 

Ji :=A[i- l/2,i + l/2) 

Jf := e~ 1/2 Ji = J v ~]^j centered at ji := iAe~ 1/2 , i e Z. 

The following lemma takes care of the case when the have not still crossed 
each other. 
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Lemma 2.2. For integers < • • • < i(n), 

(2.8) lim P(y* GJ? (fc) ,fc=l,...,n)= II P (** G J *(*))' 

where (A/,)/, are as in Theorem 2.1. 



Proof: It follows from the definition of Y t k and (2.3). □ 

Now we study the case of unordered Xj? . We consider first the case n = 2. Take 
arbitrary integers i(l) > z(2). For < 5 < 1/2 let 

Jf ■= (i ~ \ + i + \ ~ s) Ae" 1/2 = (j. + SAe~^, Ji - 5Ae~^) . 
The interval Jf ,s is strictly contained in Jf. Let 

= (^i(i(l),i(2))-l-(5,Vi(<(l),<(2)) + l + (5) A £ " 1/2 

where 

(2.9) Vi(i(l),i(2)) = ^(*(l),i(2)) = z(l)^^ + z(2)^^i 

This corresponds to the coordinates of the vector defined in (1.18) because i(l) > 
i(2). 

Lemma 2.3. For integers > z(2), and for all fixed A, 
(2.10) 

lim lim P (x k T G Jf/* , fc = 1, 2) = lim lim P (x* G Jf/* , i; fc G 7 s /* fc = 1, 2) 

5^0 e^O V T£ *( fc ) / 5^0 e^O V Te l ( fc ) Te *( 1 ).H2)' ' J 

Proof: For 7, J C K we say 

I < J if x < y for all x G 7 and y G J. 
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Considering, as before, r £ = r e — as x / 2 , let us define the intervals 

Jf(k),a = ilt(ky^i(k)) = 6-ae- 1 / 2 (l-p k _ 1 - Pk ) + x £ Ji(k)i 

J m, a = tfw + ste-wzw-ste- 1 '*) = (?_ oe -v a(1 _ Pfc _ 1 _ Pfc)+jB .^ ) 

which we get at time r £ by translating 6 Xe Jf^ and 9 Xe J e ^ at velocity — (1 — Pk-i ■ 
Pk) backwards from time r e to time r £ . Let a big enough such that 



e 

Of - 



From the one-shock law of large numbers for Xt and the definition of Y t h (which 
coincide with the for t < ti in that defintion), we have: 

Hm P(y*. G J? w>0 ,fe= 1,2) = limPpC. G J? (fc)>a , fc = 1, 2) 

= ]imP(X* G J- (fc) ,fc=l,2) 

This gives the localization of the F particles at time r £ . The idea is that those 
particles will follow their respective characteristics and thus meet during the time 
interval (r", r e ), where they coalesce and follow the new characteristic, thus ending 

in x s + 7f^ ^ at time r e . Since we are in the same scale for time and space, the 
result will follow from the law of large numbers. 

To make this rigorous, at time r" put a {a 1 discrepancies (second class particles) 
at +5Ae _1 / 2 ] + 1 and [jf^ — 5Ae _1 / 2 ] and a l \a 2 discrepancies (third class 

particles) at [7^ + tfAe -1 / 2 ] + 1 and \jf^ - 5Ae -1 / 2 ]. Label their positions Fj, 
Fj, F 2 and F 2 , respectively, for £ > r". 

Let Tg = r" + 5e -1 / 2 , where r" is defined by £F*„ = EY T n. At time r^, which 
will be smaller than r e for e sufficiently small, put a°\a 2 discrepancies at 

[j? w + (l-p - pl )(T' £ -T?)} and [^ (2) + (l- Pl -p 2 )(^-T«)]. 
Label their positions F t and Y t , respectively, for t > r' £ . 
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Since p = 1, we easily check the facts that for t > r", 

{y% e fc = i. 2 } = {*J? e = i, 2} c {yj < < F t fc } 

which shows in particular that for all t>r' e : 

(2-H) {y t \ g ^ )>a , = 1,2} c {x t 2 < ^} u {£, < F*,} 

and 

(2.12) {X* < X}} c {x, 2 < y* < X, 1 , fc = 1,2} . 

Now, a simple geometric argument relying on the laws of large numbers for 
Y^Y^Xj* (cf. Theorem 3.1 in [FF] and Remark following eq. (2.6)), (2.11) and 
(2.12) proves that 
(2.13) 

lim P (y T * e J? ( ° )ta , k=l,2) = lim P (Xi„ E J# )>a , Y T , < Y* <Y T ,,k=l, 2) 

plus 0(5), as soon as 5 > is small enough. 0(5) comes from the probability of 
the event on the right of the union on the right hand side of eq. (2.11) and the law 

of large numbers for Y_] and Y t . 

The result now follows from an application of the laws of large numbers for Y_ t , 
Y t and the fact Y t <Y t < Y t for t > r' £ . (See figure 3.1.) □ 

(Figure 3.1 enters here) 

The same argument with a more complicated notation shows the above extends 
to the case of general n, summarized below. 

Lemma 2.4. For distinct integers i(l), . . . , i(n), and for all fixed A, 



(2.14) 



&o^o P (^- Gj !(*)' fc = 1 '-"' n ) 

= J™ H P e J m ' e I -it) ,k = l,...,n), 
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where j(k) = 4>k(i(l), ... ,i(n)). 

Proof of (2.4). Follows from Lemma 2.4 and (2.3). □ 
§3. Proof of Theorems 1.1 and 1.2. 

Proof of Theorem 1.2. Let us start by noticing that Theorem 1.2 was stated in 
terms of the evolution rjt, which is not exactly the rj' t we are considering, for which 
initially the sites [c^e -1 ], k = 1, . . . , n, are occupied by second class particles. On 
the other side, it is immediate to see that coupling r\t and rj' t in the usual way, 
one concludes that the eventual discrepancies, at most n, behave as second class 
particles. They might even annihilate one another, but in any case diffuse as e -1 / 2 . 
This shows that Theorem 1.2 will be proven (i.e. for i] t ) once we check the analogous 
for our perturbed process rj' t . For this, it is sufficient, according to Daley and 
Vere- Jones (1988), Proposition 9.1. VII, to show that for every bounded continuous 
function $ with compact support J $ dA £ converges weakly to J $ dA as e — > 0. 
The former integral is 

(3-1) e 1/2 £/0Wl<)*(e 1/2 *)- 

Let M > be an integer such that the cylinder function / depends only on the 
coordinates in {— M, . . . , M}. Then, we use (2.2) to condition on the (standardized) 
locations of the shocks at time r e , and decompose (3.1) as 
(3-2) 

n 

fc=0 x6Z 

where Y £ k = e 1 ' 2 (Y^ — [x e ]) and the random variable Z s satisfies 

n 

E(\Z £ \) < H/IU J2 £l/2 H ^ 1/2 x)p(y £ fe G (e 1 / 2 (x-M),e 1 / 2 (x + M)) . 

k=0 x6Z 

By Theorem 2.1, (Y" 1 , . . . , YJ 1 ) converges weakly to ij>(X\, ■ ■ ■ , X n ) as e — > 0, and 
this implies at once that lim e ^ E(\Z e \) =0. 
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The result follows from this. To see it, notice that, for e > 0, the expression in 
(3.2) is a function of (o^J and (Y £ k ), which we denote F e {a\^ . . . , a™ , , . . . , YJ 1 ). 
By the product structure of tr^ and the law of large numbers, for all (yi, . . . , y n ) 

F £ (al e ,...,(j? e ,y 1 ,...,y n ) —> F( yi ,...,y n ) := J $ dA( yit _ tVn) 

almost surely, where ^( yi ,...,y n ) is defined as in (1.24). Both F and ip are continuous 
in (y 1; . . . , y n ). Also, {F e , e > 0} is equicontinuous as a family of functions of 
(yi, . . . , y n ). This almost sure convergence and continuity properties of F e , F and 

V> plus the weak convergence of (Y £ ) yield the weak convergence claimed in the 
statement of the theorem by standard arguments. We leave the details to the 
reader. □ 

Proof of Theorem 1.1. It suffices to check that 

n 

lim£/(0 [r . e -i +ae -i/ 2] 7k. e -i) ^ ^y Pk f p{y k <r <y k+1 ) 

where 3^fc are as in the statement of Theorem 1.1, r\t represents the process starting 
with |U e and / is a cylinder function, increasing (for the usual coordinatewise order). 
Now, to recover the above expression out of Theorem 1.2 we may take the expec- 
tation E J $<iA e where the test function $ defined by $>(w) = ^l{a < w < a + u} : 
with u > will give an upper bound for E f (6*[ X£ + a e- 1 / 2 ] r ?r e ) • Similarly by using $ 
defined by $(«;) = ^l{a — u < w < a}, with u > 0, we get a lower bound. (The fact 
that these particular test functions have two points of discontinuity is not a prob- 
lem, due to the continuity of limit measure A.) Here we are using the attractiveness 
of the system and the monotonicity of the initial profile to get fi s S(t)f < Q\[i & S(t)f ', 
for increasing continuous /. Letting u tend to zero both terms will converge to the 
desired expression. □ 

Remark 3.5 . A similar analysis can be employed to determine that the (micro- 
scopic) measure seen from ([x e + ae~ 1 ^ 2 ],r £ + se~ 1 ^ 2 ), with a, s G R fixed, converges 
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as e — > to a mixture of ^ Pfe . For that, as in the definition of ip, let 

_ f + S ( 1 ~ Pk-l ~ Pk) l{s<-t(x) 

k ~ I bk(x,t(x) + s) if s > — £(x) 

This is the function that enters in the corresponding statement. 

Remark 3.6 . In this paper we have treated in detail only the totally asymmetric 
case, where p = 1. It is not hard, but technically more cumbersome, to extend the 
analysis tol/2<p<l. One can define n' t in the same way. The extra difficulty 
comes from the fact that now the shocks can uncross each other. However, since 
the gaussian fluctuations and the law of large numbers remain valid, essentially the 
same analysis applies, with the appropriate change on parameters. Details are left 
to the reader. 
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